Electronic Structure of ZnCNi 3 
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According to a recent report by Park et al, ZnCNi3 is isostructural and isovalent to the super- 
conducting (T c ~ 8 K) anti-perovskite, MgCNi3, but shows no indication of a superconducting 
transition down to 2K. A comparison of calculated electronic structures shows that the main fea- 
tures of MgCNi3, particularly the van Hove singularity near the Fermi energy, are preserved in 
ZnCNi3. Thus the reported lack of superconductivity in ZnCNi3 is not explainable in terms of T c 
being driven to a very low value by a small Fermi level density of states. We propose that the lack 
of superconductivity, the small value of the linear specific heat coefficient 7 and the discrepancy be- 
tween theoretical and experimental lattice constants can all be explained if the material is assumed 
to be a C-deficient a-ZrCNi3 similar to the analogous non-superconducting phase of MgCNi3. 



I. INTRODUCTION 

The appearance of superconductivity near 8 K 
in the Ni-rich perovskite, MgCNi3, has stimulated 
much interest not only because it is unusual in a 
compound that is primarily Ni, but because the ex- 
act nature of the superconducting state and its mi- 
croscopic origins are still being debated. Like the 
other unusual new superconductor, MgB2, it has so 
far resisted efforts to increase the critical temper- 
ature significantly by chemical substitution. Both 
Cu and Co doping on the Ni site reduce the critical 
temperature (T c ), predictably due to band effects 
(electron doping) in the former case and possibly 
due to spin fluctuations in the lattes 2 ^. The transi- 
tion temperature can be raised by f K through Ni- 
site doping^, with Fe, but this temperature occurs in 
MgCNia-^Fe^ with x = 0.05, and any further doping 
again reduces T c . Mg deficiencies or excesses have 
some effect on the sharpness and onset of the super- 
conducting transition, but the optimal composition 
still results^' 5 in a maximum T c of 8 K. The super- 
conductivity of MgCNi 3 seems most sensitive to the 
carbon site occupancy. Boron doping on the carbon 
site^ reduces T c for relative B/C concentrations of 
up to 0.07 and eliminates superconductivity for any 
greater concentration. MgC^Nia with x <1.0 re- 
mains a cubic perovskite but undergoes an isostruc- 
tural transition 1,5 to a smaller volume a-phase that 
no longer superconducts. 

Recently the synthesis of ZnCNi3 has been re- 
ported by Park et a&. Since ZnCNi3 is very sim- 
ilar to MgCNi 3 structurally, and (as we will show) 
electronically, the lack of a superconducting transi- 
tion down to 2 K is quite unexpected. Understand- 
ing why superconductivity is seen in one compound 
but not the other could be important in resolving 
remaining questions about the unusual behavior of 
MgCNi3. The experimental data suggests^ that a 
strongly depressed density of states (DOS) (com- 
pared to MgCNi 3 ) at the Fermi level (E_p) could be 
responsible for pushing the transition temperature of 
ZnCNi 3 below 2 K. The results of a careful compari- 



son of the electronic structure of the two compounds 
are presented here, and the required lowering of the 
DOS is shown to be absent. Because the reported 
lattice constants differ by 4%, we calculate the theo- 
retical equilibrium lattice constants and explore the 
effects of pressure on the electronic structure of each 
material. We find that discrepancies when compared 
to MgCNi 3 in DOS, in lattice constant, and in ob- 
servations of superconductivity can be understood 
if the reported ZnCNi3 samples are C-deficient a- 
phase as are MgCNi 3 samples with C deficiency. 



II. CALCULATIONAL METHODS AND 
ELECTRONIC STRUCTURE 



ZnCNi3 has the typical ABO3 cubic perovskite 
structure, but with the oxygen atoms on the faces 
replaced by Ni atoms. As Zn and Mg both have a 
formal valency of 2 + , Z11CN13 is isovalent as well as 
isostructural with MgCNi3, both residing in space 
group 22f (Pm3m). Calculations were carried out 
using Wien2k 7 , a full-potential, augmented plane 
wave + local-orbital method, and with the local den- 
sity approximation (LDA) of Perdew and Wang£ to 
the exchange-correlation potential. The density is 
well-converged with 816 k-pts in the irreducible Bril- 
louin zone. The sphere radii used were 2.1 a.u. for 
Zn/Mg, and 1.72 a.u. for both C and Ni; the Rkmax 
was set to 7.00. The experimental lattice constants 
were used in the initial calculations for both MgCNi3 
(a = 3.81 A) and ZnCNi 3 (a = 3.66 A). Compres- 
sion and expansion percentages are given in terms of 
these experimental values. 

The electronic structure of MgCNi3 has been pre- 
sented previously by several group qW 1 ' 12 . The 
dominant feature is a remarkable, sharp van Hove 
singularity 65 meV below E^?, which was traced to an 
extremely flat band around the M point (=(1,1,0)—) 
of the Brillouin zone. The electronic structure of 
ZnCNi3 is very similar to that of MgCNi3 . The sharp 
peak just below the Fermi energy is still dominant, 
though it is shifted slightly downward in energy by 
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A. Pressure Dependencies and Bulk Modulus 

Expansion of the Z11CM3 lattice narrows the peak 
and brings it nearer the Fermi level (Fig. 1 middle 
panel) that is, expansion makes it more MgCNi3- 
like. This has the effect of raising the DOS at the 
Fermi energy, but the change is very small even for 
fairly large expansions. In ZnCNi3, an expansion of 
~ 12% by volume caused a change in the DOS of 
only 10%. MgCNi3 seems to be even slightly less 
sensitive than this. A calculation of ZnCNi3 at the 
equilibrium lattice constant of MgCNi3 shows, how- 
ever, that the differences in the electronic structures 
of the two compounds are due to more than sim- 
ply volume. (See bottom panel of Fig. 1). The 
probable source of the small differences in electronic 
structure is residual hybridization of Ni-d and Zn-p 
orbitals. As the lattice is expanded the overlap be- 
tween these orbitals decreases, but does not disap- 
pear completely. By looking at the orbitally-resolved 
character of the flat band, we have observed that 
the amount of Zn-p character in ZnCNi3, although 
minor, is larger by nearly a factor of two than the 
amount of Mg-p character in MgCNi3. 



FIG. 1: Top panel: ZnCNi 3 and MgCNi 3 at their 
reported equilibrium lattice constants. Middle panel: 
ZnCNi3 at its own lattice constant and at that of 
MgCNi3. The effect of pressure is rather small. Bot- 
tom panel: MgCNi3 and ZnCNi3 both at the same lat- 
tice constant (that of MgCNi3). Differences in electronic 
structure are larger than can be accounted for by pres- 
sure alone, though still small. 



approximately 30 meV and has broadened some- 
what (Fig. 1, top panel). In both compounds, the 
Ni ions are two-fold coordinated with their nearest- 
neighbors, the co-planar C ions. Hybridization be- 
tween Mg/Zn and Ni ions is very small, consis- 
tent with the very similar electronic structures of 
the two compounds. The dispersion created by the 
two-dimensional bonding of Ni-d and C-p orbitals 
is responsibl o 9 i 10 for the nearly dispersionless band 
centered on M. In ZnCNi3, the situation is much the 
same, but the 4p states of the Zn ions, with which 
the Ni ions are four-fold coordinated, do participate 
weakly in the bonding states near the Fermi energy. 
This weak bonding is three-dimensional, accounting 
for the slightly increased dispersion around the M- 
point as well as the lowered energy of the DOS peak. 

The downward shift of the peak has the effect of 
reducing the DOS at the Fermi energy (N(0)), with 
respect to that of MgCNi3, by about 1 eV -1 , i.e. 
by about 20 %. This decrease relative to the Mg 
compound is much less than what is necessary to 
account for the lack of superconductivity through 
conventional BCS theory (see Section C). 
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FIG. 2: The bulk moduli and equilibrium volumes of 
MgCNi3 and ZnCNi3. The experimental volumes are 
shown as vertical lines - that of ZnCNi3 is smaller than 
the calculated value. 

The calculated equilibrium volume of ZnCNi3 
compares very favorably with the reported value, the 
latter being 0.53% smaller by lattice constant. How- 
ever, this result is somewhat unusual in that the the- 
oretical value is actually larger than the experimen- 
tal one. The calculated equilibrium value of MgCNi3 
is 2% smaller in lattice constant than the experimen- 
tal value. The discrepancy is in the more common 
direction, but larger than usual. The two energy vs. 
volume curves are shown in Fig. 2, along with the 
experimental volumes. The bulk modulus of ZnCNi3 
taken at the theoretical volume is 251 GPa whereas 



that of MgCNi3 taken at its own theoretical volume 
is 214 GPa. The MgCNi 3 result is similar to that 
obtained through LMTO calculations^. Both these 
values are obtained by fitting the Murnaghan equa- 
tion of stated to an energy vs. volume curve and 
extracting the bulk modulus through the relation: B 

= V^|. The smaller volume Zn compound is, as 
expected, harder than the Mg compound. However, 
the bulk moduli of these two compounds, calculated 
at their respective experimental volumes, differ by 
only 3.6%, with MgCNi 3 being harder than ZnCNi 3 . 
This is a consequence of finding the theoretical equi- 
librium value of ZnCNi 3 above the reported value, 
while that of MgCNi 3 is below the reported value. 



B. Electron-Phonon Coupling 

The empirical value quoted for the electron- 
phonon coupling constant, A, in MgCNi 3 depends 
on the method by which it is obtained. Using the 
size of the jump in specific heat at the superconduct- 
ing transition and assuming weak coupling BCS be- 
havior yields^ A = 0.79. This method is obviously 
unavailable for ZnCNi 3 since no transition has yet 
been observed. A can alternatively be derived by 
the more common method of comparing experimen- 
tal and theoretical results for the linear coefficient 7 
from specific heat data: 

A = 2^-1; ^ = ^^(0) (1) 
7th 3 

There is some variation in the reported values of 
7exp for MgCNi 3 . Some sourcesiSii^ place the value 
at around 29 mJ/mol K 2 for the zero-field value, 
while othersi^i^ cite a higher value of about 33.5 
mJ/mol K 2 . Using this range of values, we obtain 
A = 1.5-1.75 for the Mg compound, in agreement 
with previous results derived in this way^ii How- 
ever, using this methodology, a negative A results for 
ZnCNi 3 , due to the small 7ex P =d-77 reported^, less 
than 25% of that of MgCNi 3 . This unphysical re- 
sult highlights the discrepancy between experimen- 
tal and theoretical comparisons of these two com- 
pounds. Large differences in observed specific heat 
data combined with very small differences in cal- 
culated electronic structure properties produce this 
unphysical value for A. 

Furthermore, it is the ratio of the 7's from the two 
different compounds that stipulates that the DOS of 
the ZnCNi 3 sample must be significantly lower than 
that of MgCNi 3 . Park et afe use the definition in Eq. 
1 along with their specific heat data to put an upper 
bound on the value of the DOS of the Zn compound 
relative to that of the Mg compound at the Fermi 
level. 



TABLE I: Comparison of MgCNi 3 and ZnCNi 3 (Exper- 
imental values taken from Park- et at) 
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It is clear from this equation that N(0)z„ takes its 
greatest value when \z n = 0. The ratio then yields 
N(0)zn < 0.41N(0)m s . As mentioned above, there 
is no such large depression of the ZnCNi 3 DOS as 
compared to the MgCNi 3 DOS at the Fermi level. 
In fact, our calculated value of N(0)z„ exceeds the 
derived upper bound by almost a factor of two. 



III. DISCUSSION 

MgCNi 3 has now been studied fairly extensively, 
and we review some results that may be relevant. 
According to Ren et al£, the carbon occupancy of 
MgCNi 3 is sensitive to preparation conditions and 
two different phases of the compound emerge. The 
a-phase is carbon depleted, while the (3 phase is 
nearly stoichiometric (carbon occupancy is 0.96). 
Both a and (3 phases share the same cubic space 
group, but the a phase lattice parameter is 1.3% 
smaller and unlike the /3-phase, it does not super- 
conduct. This is consistent with previous studies^ 
which found that T c decreases linearly with decreas- 
ing carbon concentration until eventually, at a car- 
bon occupancy of around 0.88-0.89, a multi-phase 
region is reached in which bulk superconductivity 
no longer exists. The reported a-phase occured at 
an occupancy of 0.75 at the carbon site^, well within 
this multi-phase region. Shan et a.' 16 found that the 
specific heat 7 was 50% lower in the a-phase than in 
the /J-phase. The a-phase can then be distinguished 
from the /3-phase in three important aspects: it does 
not superconduct, it has a significantly smaller 7, 
and it's equilibrium lattice constant is 1.3% smaller. 

Most if not all of the evidence regarding ZnCNi 3 
can be reconciled if we suppose that the phase re- 
ported by Park et al is a carbon deficient "a - 
ZnCNi 3 " phase corresponding to a-MgCNi 3 . The 
electronic structure of the two compounds is so sim- 
ilar that it is reasonable to assume that carbon de- 
ficiencies in the Zn structure would have much the 
same effect as carbon deficiencies in the Mg struc- 
ture. Assuming that the experimental results for this 
compound were taken from an a-phase of ZnCNi 3 , 



all discrepancies between theory and experiment dis- 
cussed in this paper disappear. A 1.3% increase in 
the lattice parameter would result in the common 
situation in which the theoretical value is smaller 
than the experimental one. If 7 is multiplied by a 
factor of two, as it would be in moving from an a 
to (3 phase, the A value calculated using Eq. 1 has 
a value of 0.42, eliminating the non-physical neg- 
ative result. Previous electronic structure calcula- 
tions show that in MgCNi3, N(0) decreases dramat- 
ically as C concentration decreases^, resulting in 
suppression of the superconducting transition. Sim- 
ilar effects would be expected for ZnCNi 3 . 

Another striking difference in the data between 
the Zn and Mg based compounds is a sharply in- 
creased lattice stiffness 0£>. Even in the minimally 
doped alloy Mgo.85Zn .i5CNi 3 , an increase of 38% 
was observed^ for 0£> and in the fully Zn substi- 
tuted compound, the increase is 67%£. The addition 
of Zn in any concentration causes a volume contrac- 
tion and a concurrent hardening of phonon modes 
in general. In pure MgCNi3, the frequency of a very 
soft acoustic Ni-based phonon mode is calculated 
in the harmonic approximation to become negative 
along much of the T — M direction of the BZ 20,21 . 
Anharmonic stabilization of this mode results in ob- 
served dynamic displacements of the Ni ions perpen- 
dicular to the Ni-C direction^. This "breathing" 
distortion allows each Ni ion to move away from its 
two C neighbors and toward the empty interstitial 
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site. In a C deficient compound, stress on the Ni 
ions would be partially relieved by vacancies, reduc- 
ing the advantage of such distortions and thereby 
increasing <->£>. 



IV. CONCLUSIONS 

From our calculations, stoichiometric ZnCNi3 and 
MgCN3 are very much alike in both structural and 
electronic properties. The experimental report of 
widely differing specific heat data and the lack of su- 
perconductivity down to 2 K seems highly unusual 
in light of the close similarity of these two com- 
pounds. The rather large suppression of the DOS 
at the Fermi energy required to interpret the exper- 
imental results using BCS theory fails to materialize 
from the calculations. All results are in line with a 
ZnCNi3 phase that is carbon-deficient rather than 
stoichiometric. Carbon deficient MgCNi3 is known 
to have a smaller volume than the stoichiometric 
compound, to have a strongly depressed 7, and to be 
non-superconducting. Our results suggest that the 
lattice constant of stoichiometric ZnCNi3 is likely to 
be larger than that which is reported (probably near 
3.74 A), and that a depression of less than 20% in 
N(0) occurs. A truly stoichiometric ZnCNi3 com- 
pound would likely be superconducting at only a 
somewhat lower temperature than MgCNi3. 
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